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This paper discusses the systematic design of an adaptive feedback linearizing neu- 
rocontroller for a high-order model of the synchronous machine/infinite bus power 
■ system. The power system is first modelled as an input-output nonlinear discrete-time 

system approximated by two neural networks. The approach allows a simple linear 
pole-placement controller (which is itself not a neural network) to be designed. The 
control law is specified such that the controller adaptively calculates an appropriate 
feedback linearizing control law at each sampling instant by utilizing plant parame- 
ter estimates provided by the neural system model. The control system also adapts 
itself on-line. This avoids the requirement for exact knowledge of the power system 
dynamics and full state measurement as well as other difficulties associated with im- 
plementing analytical input-output feedback linearizing control for a complex power 
system model. Furthermore, a departure is made from the 'ad hoc' manner in which 
many neural controllers have been designed for power systems; the approach used here 
has foundations in control theoretic concepts of adaptive feedback linearization and 
pole-placement control design. 

Simulation results demonstrate the performance of this controller for a representa- 
tive example of a single-machine/infinite bus power system configuration under various 
operational conditions. 



*This work was supported by the Natural Sciences and Engineering Research Council of Canada via a 
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1 Introduction 

The differential geometric technique of state feedback linearization has been explored over 
the past two decades for the purpose of power system control (see, for instance, [0], |L7 



IH, @])- The main idea is to perform a co-ordinate transformation of the nonlinear state 
space system equations and define a new control input so that in the new co-ordinates, the 
nonlinearities in the plant are either wholly or partially masked. This formulation may result 
in linearizations which are valid for large practical operating points of the system, as opposed 
to a local Jacobian linearization about an operating point. 

Feedback linearization has been applied to the control of power systems represented 
by the state-space single machine/infinite bus model in two main ways. The input-state 
feedback linearization approach is formulated such that the system state becomes a linear 
function of a new control input and a new state, while the output (for our purposes, the 



terminal voltage) is still a nonlinear function, e.g. see ||17|| . Although this works well for 
designing stabilizing controllers, the nonlinearities in the output map make it very difficult 
to achieve good tracking of reference voltage signals. An alternative to this approach is the 
input-output feedback linearization in which the terminal voltage becomes a linear function 
of a new control input, e.g. see f2|, ||, Due to the mathematical complexity of the 
nonlinear state-space model representing the high-order single machine/infinite bus power 
system (a 7th order model), a simplified 3rd order model is used in |2j, ||. An excitation 
controller/stabilizer is then designed and tested using the high order model of the plant. The 
effects of the unmodelled system dynamics and possible shifts in parameters of the plant raise 
some questions on the effectiveness of the controller. An inherent drawback of the feedback 
linearization approach is the non-robustness due to the need for exact knowledge, both in 
terms of the structure and parameters, of the model of the power system. Adaptive control 
has been used to compensate for parameter variation in a framework that allows the controller 



to learn the nonlinearities on-line 12 . However, this assumes that the nonlinearities can be 



parameterized linearly in some unknown parameters ||12|| , [|T3|]. Some power system stabilizers 
based on adaptive control design techniques have been shown to increase the operating range 
over which they can provide good control as well as exhibit robustness to parameterized 
system disturbances . 



This paper presents a design which addresses these problems by working directly with the 
more accurate 7th order model of the power system. The nonlinear dynamical system is mod- 
elled using multilayer neural networks made of sigmoid type nonlinearities and adjustable 
weights which are nonlinearly parameterized. The flexibility provided by the nonlinear pa- 
rameterization allows for a more realistic representation of the underlying nonlinear power 
system model. Power system controller designs using neural networks of various configura- 
tions exist in the literature and are reported in works like p|. Neural adaptive control of 
feedback linearizable nonlinear systems was first proposed in and extensively analyzed in 
||. In these references, the networks were trained using the backpropagation algorithm ||; 
the nonlinear function approximation capability of neural networks (see jnj) was exploited 
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to model the plant in order to determine an appropriate feedback linearizing control. How- 
ever, the approach used in J/J, || provides no clear way of specifying closed loop dynamics. 
In this work the neural network identifier is trained in batch form and adapted using the 
Levenberg-Marquardt optimization (as in |15| ) rather than using standard backpropagation. 
Also, the whole control problem is formulated such that tracking dynamics can be user- 
specified by approprately assigning closed loop system poles. Weight adaptation based on 
nonlinear optimization is attractive here because it significantly corrects some of the defects 
of standard backpropagation, most importantly slow rate of convergence. In this framework, 
the trained neural networks give an on-line estimate of the parameters of the power system 
at every sampling instant which are then used to calculate the feedback linearizing control 
law. The controller is designed to provide reference tracking (with bounded state), damping 
of power angle oscillations and tolerate variations in the plant parameters. It may also be 
adapted on-line. 

This paper is organized as follows: Section 2 presents the state-space synchronous ma- 
chine/infinte bus power system model. The system relative degree and minimum phase 
properties, important for feedback linearization to be successfully applied, are also discussed. 
The class of systems to which this model belongs and which motivates the choice of a neural 
identifier is determined. Section 3 focuses on the neural modelling of the power system while 
Section 4 presents a systematic design of the proposed excitation controller /power system 
stabilizer. Finally, simulation results are presented and discussed in Section 5, followed by 
some concluding remarks. 



2 System and Identification Model 

The model used is developed in |TjJ (p. 54ff). It yields a Park's 7th order nonlinear time- 
invariant, state-space description of the synchronous generator /infinite bus system. This is 
connected through a balanced pair of transmission lines to a power system modelled simply 
as one of infinitely large ability to supply and absorb power at a set voltage level, i.e. an 
infinite bus system. The system is depicted in Fig. [I]. 
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Figure 1: Synchronous Machine Connected to an Infinite Bus 



State variables for this generator /infinite bus set-up are selected as follows: 



• 5(xi) - power angle 

• uj{x 2 ) - power angle derivative 

• \d{%3) - d-axis flux linkage 

• A g (s 4 ) - q-axis flux linkage 

• Xf(x 5 ) - field flux linkage 

• Xk d (xQ) - damper d-axis flux linkage 

• Xk q (x-j) - damper q-axis flux linkage 

The system model is given by 

5 = us 

03 = ^[P m -Pe(S,X d ,X q )-Du] (1) 

Zn 

X = u b ([RL- 1 + Z]X + [v d v q v f 0) T ) 

where A := [Ad X q A/ Xk d Xk q ] T , are the per unit flux linkages, R := diag [r s r s — r/ — r\ d — r^J, 
are per unit resistances. H is the inertia constant, D encapsulates the damping in the system, 
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As shown in Fig. [I], the control input {u) is the field voltage Vf and the output to be 
controlled (y) is the generator terminal voltage Vt given by 



v t = h(5, A d , X q ) = \Jv 2 d + v t 



where 



)' d 
v„ 



r u I d — x u I q + v oa [Acos(n/2 — 5) + Bsin(n/2 — 5)] 
r\\l q — X\\ld — v 00 [Bcos(iy/2 — 5) — Asin(ir/2 — 5)]- 



the variables are described in |L6| and [Q. A, B are constant matrices having values de- 
pending on physical parameters of the power system, Voo is the voltage of the infinite bus, 
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are d— and q— axis currents, rn, iy are components of transmission line resistance and 
reactance (see jnj). The system model ([I]) belongs to the class of control affine nonlinear 
systems of the form 

x = f(x)+g(x)-u (2) 

y = h(x). 

Here /, g are smooth functions; the state x G M. n , input u G R and output y G R. The 
relative degree and minimum phase properties of this system are examined next. 

The relative degree of the power system ([]]) may be obtained by differentiating y (i.e. 
Vt) until u (vf) appears. This is achieved by computing the Lie derivatives (see @) of the 
system. 

Lfh = Vhx 

dh dh dh 

= -j—xi + • • • + 3—^5 + • • • + -1—X7 
ax 1 dx§ dx-j 

Observe that Vf appears explicitly in the expression for x 5 (Xf) in (jH). Therefore the output 
need only be differentiated once before the control input Vf appears, i.e. the relative degree 
is 1. It is easy to show that the system (^) is nominally minimum phase over the practical 
operating region by linearizing the system at various operating points in the state-space and 
verifying left half plane zeros (see Appendix B in 



From the relative degree and minimum phase properties, the identification model pro- 
posed for the system will be the single-input single-output relative-degree-one nonlinear 
discrete-time system 

y[k + l] = f[y(k),... ,y(k-n + l), u(k-l),... ,u(k-m)\ + (3) 
g[y(k), ■ ■ ■ , y{k - n + 1), u(k - 1), . . . , u(k - m)\ ■ u{k) 

where y is the output, u is the input, / and g are assumed smooth and the system is assumed 
minimum phase. For the 7 t/l -order power system with relative degree 1, n = 7, m = 6. 



3 Modelling the power system 

The first stage in the controller design is to identify the synchronous machine/infinite bus 
dynamical model using the structure (§). The set-up to do so is shown in Fig. ||. TDLs are 
tapped delay lines used to implement the regressor structure for both networks. NN-1 and 
NN-2 are feedforward neural networks whose input vector consists of present output and 
delayed inputs and outputs. They are used to predict the next output yik + 1). iV7V-l has 
one hidden layer made up of a total of p neurons with hyperbolic tangent tanh activation 
function and an output layer of one neuron with linear activation function. NN-2 has q 
tanh hidden layer neurons and one output neuron. All weights are initialized randomly and 
contain biases. 

A training set is generated by exciting the plant with a variety of input signals u (of ap- 
propriate range) and measuring the output y. That is, let Z N = {[u(k),y(k)]\k = 1, . . . , iV}, 
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Synch, machine model 




Figure 2: Neural system identification of the synchronous machine model 



represent the training data pairs and N be the total number of data sets. The estimate of 
the plant output given by the neural nets is 

y[k + 1] = f[z(k), w] + g[z(k),v] ■ u(k) (4) 

where 

p / m+n \ 

f(z(k),w) = ^^Wi tanh I WijZj + bi J + b 

i=l \ j=l ) 

p / m+n \ 

g(z(k),v) = tanh I ^ v i3 z 3 - + J + a; 

i=l \ j=l ) 

(w,b) and (v,a) are weights and biases respectively for the /(•) and g(-) networks. Let the 
set of possible weights be 9 and define the cost function J as the mean square error in 
predicting the plant output. That is, 

1 N 

m®, zN ) = ^J2 + !) - M + 1 i )] T + x ) - y( k + ii )] • ( 5 ) 

fc=i 

The optimal weights are determined by solving 
= arg mm Jn(Q, Z n ) 
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Let the weight update rule toward achieving the optimal choice be given by 

©fc+l = ©fe + PkSk 

where ©& is the current iterate of the set of weights, pk is the step size and Sk is the search 
direction. 

For identification purposes, the weight update is done in batch form using the Levenberg- 
Marquardt optimization as described in When the cost function is minimized to an 

acceptable level, weight adaptation is discontinued and the neural model is validated on the 
training data set and subsequently cross-validated on other sets of data not in the original 
training set. Simulation results are provided which show the effectiveness of the proposed 
identification technique. 



4 Adaptive linearizing control design 

Recall that the power system is modelled by the SISO discrete time system: 



y[k+l] = f[y(k),...,y{k-n + l),u(k-l). 

g[y(k), ... ,y(k-n + l), u(k - 1), 

and estimated by the neural model: 

y[k + l] = f[y(k), . . . , y(k — n + 1), u(k - 1). 

g[y(k), ■■■ ,y(k-n+l), u(k - 1), 



, u{k — m)} + 
, u{k — m)} ■ u(k) 

,u(k — m)] + 
, u(k — m)] ■ u(k) 



(6) 



(7) 



A simple pole-placement controller is proposed to allow the tracking dynamics to be specified 
by appropriately placing closed-loop poles at stable locations. Let the reference signal to be 
tracked be given by r(k). 

To facilitate the discussion, a state-space form of the system (|6|) is constructed by choosing 
state variables as follows: 

2n (k) = y(k-n+l) 



z\ n {k) = y{k) 
Z2i{k) = u(k — m) 

Z2m(k) = U(k - 1) 

and we subsequently define 



Zl 
Z2 



Z21 



z 2i 
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This leads to the state- 
z n (k + l) = 



:-space model: 
z\ 2 {k) 



z ln (k + l) = f[z(k)] + g[z(k)] ■ u(k) = y[k + 1] (8) 
z 2 \{k + l) = z 22 (k) 

z 2m (k + l) = u(k) 
y(k) = z ln (k) 

Now we define a new control input u and set 

u{ ® = -Arr-A-mm+m]-, (9) 

gl z ( k )\ 

observe that this results in input-output closed loop dynamics of 
y[k + l] = u(k) 

Now we address the primary control objective of tracking a constant reference voltage r = v t . 
This motivates the choice of the feedback control u(k) as 

u(k) = K ir (k) - [C p - iy (k) + C p ^ 2 y(k - 1) + • • • + C y(k -p + 1)] (10) 
with p G Z + and Ci G 1R. Hence, 

y[k + l]+C p ^y(k) + C p ^y(k-l) + --- + C y(k-p + l) = K ir (k). 
Therefore the transfer function between the reference and y is given by 

YM IK = Zi _. (id 

V nf [z] R[z) zp + C^zp- 1 + ■ • • + dz + C ' Q(z)' 1 1 

Closed loop poles for the system can be assigned by selecting constants Cj's so that the zeros 
of the characteristic polynomial Q(z) are all in the unit circle. 

Some assumptions made on the system and a discussion of the stability of this kind of 
control scheme can be found in ||18|| . Specifically, it is shown therein that the approach used 
here is a more general form of that extensively analyzed in || . 

Now we rewrite the plant in the input-output form as 

y(k+l) = f[z(k)] + g[z(k)] ■ u(k) (12) 

The plant (0) is thus modelled by the neural networks as: 

y(k + l) = f[z{k),w] + g[z(k),v]-u{k) (13) 

For the three-layer neural networks used in the identification stage, with total hidden neurons 
p and q respectively and weights/biases (u> i5 ibi G w) and (t>j, ij 6 v) , the functions /(-,•) 
and g(-, •) can be expressed as 

p / m+n \ 

f[z(k), w] = tanh I ^ w ij x j + w { J (14) 

i=i \ j=i J 
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and 

g[z(k),v] = tanh ^ VijXj + Vi (15) 



=1 




Funahashi's theorem [|nj] says that we can approximate / and g as close as we like by a 
suitably parameterized neural network. It is not presently known how to choose the number 
of neurons given a desired error tolerance. However, the typical approach is to choose a con- 
figuration and try to obtain good weights/biases by trying to solve a nonlinear optimization 
problem. If the approximation is poor, the configuration is adjusted and repeated. This will 
be accomplished via an "off-line" training phase using the Levenberg-Marquardt algorithm 
(as discussed in Section ||), to obtain initial estimates of / and g. 

Now, we focus on the on-line phase, in which the controller is applied to the system. An 
additional constraint to consider is the need to ensure that g remains nonzero due to the 
fact that the controller incorporates a division by g as shown in (|16|). 

Let w(fc) and v(/c) denote the estimates of w and v at sampling instant k. The control 
input u(k) is then given as: 

Control Law: 

-f[z(k),w(k)]+u(k) 

«(«) = -r-y- tt-^ 16 

g[z(k),v(k)\ 

with u(k) as previously defined in fllCf) . The control u(k) is applied to both the plant and 
the neural network model; the network weights are updated using the error between the 
plant output and the output predicted by the neural model. At each sampling instant k, the 
estimated plant output is given by 

f(k + l) = /[z(*),w(A0] + g[z(k)Mk)]-u(k) (17) 

The error associated with this estimate {e*(k + 1)} is defined as 

e*(k + l) = y*(k + 1) -y(k + 1) (18) 

This error (together with a deadzone adaptation) will be used in the weight update rule 
described next. A deadzone is an absolute lower limit on the adaptation error below which 
no adaptation should be done. That is to say, before the weights are adapted at each 
sampling instant, the plant output error e*(k + 1) is compared with the dead zone radius d . 
If \e*(k + 1)| < do, no adaptation is done, otherwise the adaptation is done using an error 
adaptation rule based on the deadzone D(e) given by 

r if |e| < d 

D(e) = < e - d if e > d (19) 
y e + d if e < —d 

which ensures that the constraint of nonzero g holds. In order to implement the weight 
update rule, first define the set of network weights as 

w 

= 

v 
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and the Jacobian to be computed at each time step by 



dy*(k + 1) 



0(fc) 



a/[z(*0,w] 

<9w 



w(k; 



dg[z(k),v] 
dv 



v(k) 



u(k) 



(20) 



This Jacobian can be calculated using the backpropagation algorithm || at the end of each 
iteration, since the variable y*(k + 1) is the output of the neural networks, and is available 
for use ||. 
Adaptation Rule: 



9(* + 1) = 9(*) - D [e\k + 1)] J k (21) 

[i + J k Jk) 

Power system stabilizer: The power system stabilizer adds damping to the generator rotor 
angle oscillations by modulating its excitation using auxiliary stabilizing signal(s). In order 
to accomplish this, the stabilizer should be desgined to produce a component of electrical 
torque in phase with the rotor speed deviations. For the machine model being used, the 
signal of choice for modulating the excitation is the rotor angle derivative 5. This requires 
that the linearizing control law @ be modified to introduce a torque component which 
damps out the oscillations of the rotor angle. The new control law is given by 

u(k) = U{k \~J { ' ] + K pss ■ Kk) (22) 

where K pss is a suitable weighting of the rotor angle derivative. The gain K pss for the neural 
controller is selected using a root locus criterion extensively discussed in ||. The proposed 
excitation controller/power system stabilizer control set-up is depicted in Fig. |3|. 



5 Simulation Results and Discussions 

In order to accomplish the system identification part, the power system is modelled in 
SIMULINK and sampled at 2 ms to ensure the fast electrical dynamics of the plant are 
adequately captured. Next, random input signals Vf G [—0.1 0.1] are used to excite the 
plant and the terminal voltage Vt is measured for each case. 10,000 such data pairs (vf, Vt) 
are collected and divided into two parts - one for training the neural nets and the other for 
cross-validating the resulting neural model. The neural nets have the following structure 

• /(■) and g(-) both have 5 hidden hyperbolic tangent neurons with bias and one linear 
output neuron 
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Figure 3: Adaptive linearizing control/power system stabilizer set-up 



• The weight matrices for both /(■) and g(-) have dimension 5 x 14 in the hidden layer 
and 1 x 6 in the output layer. 

• The input to each network consists of current and delayed Vf and Vt up to the most 
delayed signal in each case (13 in all) with the bias term accounting for the last column 
of the weight matrix 

The optimization routine used for offline weight adaptation is allowed to proceed for a total 
of 150 iterations at the end of which the cost function is minimized to the order of 10 -6 . 
Fig. |] (top) shows the cross-validation of the model on a portion of test data not included 
in the original training set while Fig. |] (bottom) shows the prediction error on this data 
(from this, the deadzone radius d is selected as 0.01). Figures showing the validation of 
the neural model on the training data are not included since they do not, in general, give a 
reliable indication of how well the neural model performs. Observe from the error plot that 
the error is essentially close to zero for most of the samples but is more pronounced around 
the range 3000 - 4000; even in this range, the errors are not unduly large. The maximum 
value of e*(k + 1) in this range is 0.06 resulting in a small relative error of 

|f * L 006 x 100% = 2.7% 



\vL 2.25 

Observe also that this is the maximum error and it only occurs for a short sampling range. 
We argue that adapting the weights of the network on-line will serve to compensate for the 
effects of this error on controller performance. 
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Figure 4: Validating the neural identification scheme on test data 



We begin with voltage tracking only and compare the neural excitation controller with 



both a conventional ST1A high gain excitation controller (see [0 p. 365) and a nonlinear 
analytic feedback linearizing excitation controller proposed in ||. The gain and time constant 
of the ST1A excitation controller are K e = 200 and T e = 0, so that the control input is given 
by 

v f = 200^(v ref -v t ). 

All quantities are in per unit and 
^- = 3.9056 x 10" 4 

so that the control law is given by 

Vf = 0.0781(i; re/ - v t ) (23) 

The ST1A exciter is based on a linearization at the operating point v re f = 1.1392. 
The analytic excitation controller is designed in || and has a control of the form 

v = -K(v t -v ref ), K>0 (24) 
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where v is a virtual control resulting from performing input-output feedback linearization on 
a reduced order model of the power system. The closed loop transfer function is given by 



v ref {s) 



K 



s + K 



so that the terminal voltage tracks the reference exponentially at a rate depending on the 
controller gain K. For the simulations done here 



K 
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in order to provide a transient response similar to that of the ST1A controller. 
Recall that the neural excitation control law is given by 



v f (k) = 



g(-) 



where 



u(k) = K ir (k) - [C p _! C p _ 2 ... Co] 



Vk 
Vk-i 

V(k-p+i) 



(25) 



and r(k) is a constant reference voltage. The constants C are obtained from the character- 
istic polynomial of the desired closed-loop system given by 

Q(z) = z p + C p _ 1 z p ~ 1 + : . . . , +C lZ + C 

and Ki is chosen to be 1 + Y^i=o to yield unity DC gain. The transfer function is then 
given by 



V t [z] _Y[z] 



V ref [z] R[z] z'f + Cp^zP- 1 + • • • + C 
It remains to select p appropriately so that 'nice' control action is obtained. 
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ome special cases 



p=0: Observe that setting p — 0, K\ = 1 causes (|5|) to reduce to 
u{k) = r(k) 

and this is how the reference signal is obtained in the work of Chen and Khalil H. 
This approach gives a control structure which lends itself to stability analysis and 
allows parameter convergence to be guaranteed. However, the dynamics of the system 
can not be prescribed (for instance by trying to place poles of the closed loop system). 
Furthermore, using this scheme for the excitation control/pss application shows a harsh 
demand for tracking and very high level of control action. The control set-up in this 
case is shown in Fig. [5] while simulation result is depicted in Fig. |6]. 



r(k) 



r(k)-f u(k) 



Nonlinear system 

and neural 
network structure 



n 



y(k+l) 



Figure 5: Excitation Control setup for Case p = 



Case p = 1: Setting p = 1 causes ( p5|) to reduce to 

u(k) = (l + C )r(k)-C -y k (26) 

so that the transfer function of the resulting closed loop system is 

Y[z] _ l + C 
R[z] z + C 

with a first order characteristic polynomial. This should give a 'nice' first order track- 
ing behaviour and constrain the level of control effort. To check the veracity of this 
reasoning, the controller /pss is tested with u defined as in (|26|) and Co selected to give 
similar speed of response as the ST1A conventional controller used for comparison. 
The control setup in this case is shown in Fig. ^ while the response is depicted in 
Fig. H, where it is compared to the analytic controller design based on a reduced order 
model of the machine reported in the work of Kennedy, et. al. Observe that the 
neural controller is no better than the analytic controller in this case. Furthermore, it 
is apparent that the response is not that of a first order linear system. In order to do 
a fair comparison with the analytic feedback linearizing controller proposed in ||, we 
also compare the response for p = 3 (since M used a simplified 3rd order power system 
model for controller design). In this case, (to place the closed-loop poles at the same 
locations as the analytic controller, z « 0.7), the transfer function is of the form 

Y[z] K x 0.03 

R\z\ ~ z 3 + C 2 z 2 + dz + C ~ z 3 - 2.1z 2 + 1.47,2 - 0.34 
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Figure 6: Excitation Controller for Case p = 



The response is depicted in Fig. ||. The neural controller is clearly better in this case 
since the 3rd order analytical feedback lineariztion did not decouple the high order 
modes. 

Case, p = 7: the effects due to inexact pole/zero cancellation mean that we still have 
a 7 th order system to deal with. It is proposed that a controller be designed to "force" 
the system to behave like a stable 7th order system. The control setup is shown in 
Fig. 1(| In practice, we know that exact pole/zero cancellation will not result from 



feedback linearization, hence we will select p=7 and attempt to impose the desired 7 th 
order system behaviour and suppress the effects of inexact cancellation. Let us choose 
to place the closed loop poles at z = 0.7, to mimic the behaviour of the ST1A controller 
at the operating point. The resulting characteristic polynomial is given by 



Q(z) 



[z-0.7Y 



z 7 - 4.9/ + 10.29/ 



12.005/ + 8.4035/ - 3.5295/ + 0.8235,2 



-0.0824 



Simulation results shown so far suggest that this later scheme gives better performance 
than the others considered. The main challenge at this time is how to extend the 
same type of stability analysis and convergence proofs given in || to this particular 
framework or develop alternative methods of analysis. 
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Figure 8: Excitation Controller for Case p = 1 



Now we compare the response of the proposed neurocontroller (with p — 7) the analytic 
controller of || and the ST1A conventional controller. 

The system output is initially at the nominal value of Vt = 1.1392. Then a step change of 
0.1 p.u. in reference voltage occurs at t = 1.0 second. The performance of all three controllers 



is depicted in Fig. 11. At t = 1.3 seconds, = 0.9950 for the ST1A exciter, 1.0002 for the 



analytic controller and 0.9979 for the neural controller. At t — 2.0 seconds, = 0.9952 
for the ST1A exciter, 1.0000 for the analytic controller and 0.9998 for the neural controller. 
The analytic controller initially gives a slight overshoot before settling to the final value. 
The neural controller has a steady state error of 0.02% while the ST1A exciter gives an error 
of 0.5%. The control effort is just less for the neural and ST1A excitation control than the 
analytic controller. 

Next we test the ability of the neural and analytic excitation controllers to track reference 
signals at an operating point far from the nominal, i.e. from v re f = 2.0 to v re f = 2.1. The step 



change takes place at t = 1.0 second and the response is depicted in Fig. 12. It is observed 



that similar transient response characteristics are obtained compared to the response to a 
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Figure 9: Excitation Controller for Case p = 3 
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Figure 10: Excitation Control setup for proposed approach (Case p — 7) 



perturbation about the nominal operating point as shown in Fig. [□]. Also the control effort 
Vf is less for the neural controller than for the analytic controller. Finally, observe that 
in both Fig. [ll| and Fig. 



12, the magnitude of the rotor angle oscillations for the neural 



excitation controller mimics the behaviour of the system under no control in response to a 
small perturbation. This is due to fact that feedback linearization was applied to a 7 th order 
system, as opposed to a 3 rd order, as done with analytic feedback linearization. Hence, the 
effect of excitation control is almost completely decoupled from the mechanical modes of the 
system. 

Excitation controller/power system stabilizer: In ||, the control law used for the 
power system stabilizer is given by 



-K(v t - V ref + Kt8), 
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Figure 11: Excitation control-perturbation from nominal operating point 

giving voltage tracking dynamics: 

v t = -K(v t - v ref ) - KKi5 
with gains 



K 



200, 



0.7091 



These values are determined by a root locus criterion (see ||) in order to obtain optimal 
damping of rotor angle oscillations. 

Recall that the control law for the neural adaptive controller /power system stabilizer is 
given by 

u k = m ~/ { ' ] + K pss ■ 5(k) 



where the gain K pss is given by 



v x 0.7091. 



(27) 



For the simulations which follow, v — 3. In general, increasing v has the effect of improving 
the damping on 5 for step changes in reference signal. However v can not be assigned 
arbitrarily large values. For instance, values of v much higher than 3 have the effect of making 
the system unstable. The output of the system is originally at the nominal value of v t = 
1.1392. A step change of O.lp.u. in reference voltage occurs at t = 2.0 seconds. Observe from 



Fig. [13| that the neural controller accurately tracks the reference signal (actually, tracking 
error is 0.0082 % after 2.8 seconds) and quickly damps out the rotor angle oscillations. 
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Figure 12: Excitation control for change in operating point far from nominal 



As done for the excitation controller, the controller /pss is validated for a change in 
operating point to v re f = 2.0 and a step change of 0.1 p.u. takes place at t = 2.0 seconds. 
Dynamic behaviour similar to that in Fig. |13| are obtained for both controllers. The figure 
depicting the tracking and stabilizing dynamics for this new operating point is shown in Fig. 


The neural controller was also tested for big changes in the operating points of the 
system as follows: a step increase in v re f from 1 p.u. to 2 p.u. was introduced at t — 0. 
After 2 seconds, v re f dropped to 0.5 p.u. Fig. [15] depicts the results for both the neural and 
analytic controllers. Both controllers performed well and exhibited similar tracking dynamics 
compared to the response to a small step change in reference voltage shown in Fig. [T3|. 

The tuned analytic controller described above exhibits high frequency harmonics during 
tracking due to the effects of the unmodelled dynamics in the design phase. The neural model 
does not exhibit such pronounced oscillations. Also observe that the control effort (vf in 
Figs. O and Il5|) is several orders of magnitude higher for the analytic controller than for the 



neural controller. Finally, observe that the steady state tracking error is effectively zero in 
both cases and the speed of response is comparable. Another set of tests done on the control 
system involves the inertia parameter H drifting from its nominal value of 9.5 to about 4.75 
(i.e. a 50% drift) to see how well the control system behaves when feedback linearization is 
inexact. This is depicted in Fig. |16[ The neural controller performs better than the analytic 
controller in this case since the analytic controller design was based on exact knowledge of 
all the parameters of the plant. It is observed that the tracking dynamics of the neural 
controller is essentially the same as what we have in Fig. where the parameter H = 9.5. 
For this test, the value of v in (]27f) is reduced to v = 2.6. For arbitrary variation in the 



parameters, the neural controller performance can be improved by reducing the size of the 
deadzone used in the adaptation algorithm so that weights are updated for values of do less 
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Figure 13: v t , Vf and 8 behaviour for a 0.1 p.u. step change in ref. voltage occuring at 
t=2.0 seconds 



than 0.01 in this particular case. 

The system performance under a disturbance is investigated by forcing the mechanical 
power input, P m to drop suddenly by 33.4% (from 1.6512 to 1.1); the response is shown in 
Fig. [17]. Observe that both controllers recover quickly and restore the terminal voltage Vt to 
its nominally value. The rotor angle changes to a new value but still retains its damping. In 
this case both controllers gave the same level of performance, although the tracking dynamics 
for the neural controller is just slightly better. We note that in practice this may be a rather 
drastic change in power level, but the simulation is done just to see how far the system may 
be "pushed" whilst still maintaining control action. 

Finally, the effect of on-line weight adaptation on the response of the overall control 
system is investigated by observing the rate at which the error adaptation occurs. The 
Fig. [18] shows how quickly the plant output error {y{k + 1) — y{k + 1)} is adapted for the 
parameter variation, disturbance and normal excitation controller/power system stabilizer 
tests. Observe that the errors are not too large (apart from that associated with system 
disturbance) and are adapted quickly. Therefore it is concluded that adapting the weights 
online does not unduly slow down the response of the overall control system. 



6 Concluding Remarks 

A neural adaptive feedback linearizing excitation controller/ power system stabilizer has 
been proposed for a high-order model of the synchronous machine/infinite bus power system. 
The controller provides tracking performance with less high frequency harmonics than an 
analytic controller based on a reduced-order model of the system and would seem to tolerate 
parameter drift. In the cases considered, step changes in reference voltage are tracked, with 
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Figure 14: Vt, Vf and 5 behaviour for a 0.1 p.u. step change in ref. voltage occurring at 
t=2.0 seconds for operating point far from nominal 



the neural controller requiring significantly less control effort than the analytic controller. 
The neural controller also provides damping of power angle oscillations with similar voltage 
and power angle transient characteristics over different operating points of the system. Also 
on-line controller adaptation does not slow down the overall control system response since 
output errors are reduced quickly. The advantage of this approach lies in the fact that exact 
knowledge of the power system dynamics is not absolutely necessary and the system state 
need not all be measurable. Also, the complicated mathematical mapping required to carry 
out analytic feedback linearization is avoided. The control scheme is quite straightforward 
and easy to synthesize. Thorough mathematical analyses guaranteeing convergence and/or 
stability for this type of control scheme (comparable to that reported in ||) is still to be 
done. 
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Figure 17: v t , Vf and 5 behaviour under a disturbance 
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Figure 18: Plant output errors for varying parameter (top), disturbance (middle) and exci- 
tation controller /pss (bottom) 
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